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<^ I We conjecture the existence of a relationship between frustration and the transition point at 



zero temperature of the Ising spin glasses. The relation reveals that in several Ising spin glass 
models the concentration of ferromagnetic bonds is close to the critical concentration at zero 
temperature when the output of a function about frustration is equal to unity. The function 

c/3 ■ is the derivative of the average number of frustrated plaquettes with respect to the average 

number of antiferromagnetic bonds. This relation is conjectured in the Ising spin glasses with 
binary couplings on the two-dimensional lattices, the hierarchical lattices, and the three-body 
. Ising spin glasses with binary couplings on the two-dimensional lattices. In addition, the 

C ■ same argument in the Sherrington-Kirkpatrick model yields a point that is identical to the 



replica-symmetric solution of the transition point at zero temperature. 

KEYWORDS: spin glasses, phase diagrams, frustration, ground states, finite dimensions, Slierrington- 
Kirkpatrick model 



1. Introduction 

Phase transition is one of the fundamental issues in physics. Most of the successes on the 
topics in the field of spin glasses''^ have been limited to the mean-field models.^ Unveiling 
the structure of phase diagrams in finite-dimensional spin glasses remains one of the most 



H ■ challenging problems. In particular, although the phase transition to the spin-glass phase is 

a characteristic phenomenon in spin glasses, we have not had a clear understanding of its 
nature in finite dimensions. Moreover, the determination of the order-disorder phase bound- 
aries in spin glasses has a practical significance beyond a pure physicists' interest, since the 
boundaries in spin glasses correspond to the accuracy threshold in topological quantum error- 
correcting codes."* 

Concerning the phase diagram, there has been a proposition based on the property of the 
Nishimori line (NL),^'^ which is a special line across the phase diagram as depicted in Fig. 1, 
that the phase transition between the ferromagnetic and non -ferromagnetic phases at lower 
temperature than the NL is induced by a geometric nature. ^'^The form of the entropy of frus- 
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\ 
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Fig. 1. The phase diagram of the ±J Ising model in two dimensions. The vertical and horizontal axes express 
temperature and the concentration of ferromagnetic bonds, respectively. The symbols P, F, and SG denote the 
paramagnetic, ferromagnetic, and spin-glass phases, respectively. The Nishimori line (NL) is drawn as a dashed 
line. The critical concentration at zero temperature is denoted by pc- 



tration distribution is identical with the form of the free energy on the NL. Frustration is a 
geometric quantity defined later. Since phase transitions are usually signaled by singularities 
in the free energy, the entropy of frustration distribution also has a singularity at the transition 
point on the NL. This is the origin of the geometry-induced phase transition. However, this 
proposition is a conjecture since the singularities in the free energy and the entropy of frus- 
tration distribution do not necessarily coincide in other regions than the NL. The case on the 
NL, nevertheless, suggests the potential of frustration in the phase transitions in spin glasses. 
In order to illustrate the concept of frustration, we consider the ±7 Ising model on the 
square lattice governed by the Hamiltonian, 

H = -J]j,jS,Sj, (1) 

where 7,y takes either 7(> 0) with probability p or -7 with probability I - p, and 5, is 
an Ising variable (5, = ±1). The summation runs over nearest neighbors. If the number of 
negative (antiferromagnetic) coupling constants is odd in a loop composed of bonds, there is 
no spin configuration permitting all bonds to be the lower energy state {-JtjS tS j = -J). The 
product of coupling constants n<r7>ef Jij over an arbitrary loop c is called frustration.^-^ If the 
frustration of a loop has a negative value, the loop is regarded as a "frustrated loop". 

Motivated by the argument about the geometry-induced phase transition, we aim to detect 
phase transitions only with consideration for frustration in the present paper. More specifi- 
cally, we concentrate our attention on the number of frustrated plaquettes in the Ising spin 
glass models. Note that a "plaquette" means an elementary loop composed of bonds on the 
lattice, which cannot be divided into any other loops. We observe the change of the aver- 
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age number of frustrated plaquettes as the average number of antiferromagnetic bonds is 
increased. 

The present paper is organized as follows. In Sect. 2 we consider the ±7 Ising model on 
the two-dimensional lattices, and an interesting relationship between the average number of 
frustrated plaquettes and the transition point at zero temperature is found. In Sect. 3 the cases 
of the hierarchical lattices are examined. We confirm the same relationship in the model with 
three-body interactions in Sect. 4. The Sherrington-Kirkpatrick (SK) model is investigated in 
Sect. 5. The final section is devoted to summary and discussion. 

2. Ising spin glasses witli binary couplings on the two dimensional lattices 

Let us consider the frustration in the +/ Ising model (Eq. (1)) on the square lattice. We 
investigate the frustration of plaquettes, which are the squares composed of four bonds in 
the case of the square lattice. In order to understand the frustration in the system, we depict 
Fig. 2. Dashed lines express the square lattice and antiferromagnetic bonds are traversed 
by thick segments. A circle in a plaquette means that the plaquette is frustrated. The figure 
shows that ends of a segment identify frustrated plaquettes. Consider a loop on the lattice and 
a segment i. We represent the number of intersections of the loop and the segment as a, and 
the number of ends of the segment inside the loop as bj. If one of the ends lies inside the loop 
and the other end lies outside, the segment has to cross the loop odd times, that is, a, is odd. 
In contrast, if both the ends lie on a same side of the loop, a, is even. A relation 

(-ir = (-1)"' (2) 
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Fig. 2. Antiferromagnetic bonds and frustrated plaquettes on the square lattice. Thick segments traverse anti- 
ferromagnetic bonds, and circles denote frustrated plaquettes. Two segments cross a loop at the lower left-hand 
corner One of the segments crosses the loop two times. An only circle lies inside the loop. 
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is satisfied as a result. Since tlie similar equations are derived for any segment, we have 

(-1)^'"' = (-1)^'^', (3) 

where the summation runs over segments on the lattice. The left-hand side of this equation is 
the product n<i7>ec(-^(7/-^)' where c denotes the loop. Thus, we extract only the parity of the 
number of circles inside the loop from the product. When the loop contains even circles, the 
product misses them, but no spin configuration realizes the lower energy state in all couplings. 
Hence, it is appropriate to consider elementary loops, namely plaquettes, to investigate the 
frustration in the system. 

We calculate the configurational average over the distribution of {7,^} of the number N^_^ 
of frustrated plaquettes. ^'^' ^^ Note that the concentration of ferromagnetic bonds is denoted by 
p and the periodic boundary condition is assumed. We have 



K\p) ■= Wl,, = 



c ^ <ij>ec ' 



(4) 



2 
where A^ is the number of spins and [ ]p denotes the configurational average with p. The 

number of plaquettes is equal to the number of spins on the square lattice. The summation 

has run over all plaquettes c. 

We next differentiate N^^^ip) with respect to the average number Nl^{p) of antiferromagnet 

bonds, which is equal to 2A'^(1 - p) in the case of the square lattice, to calculate the change of 

the number of frustrated plaquettes by the increase of antiferromagnetic bonds. The resulting 

derivative is 

v^»:=P^ -^=2(2p-lf. (5) 

\ dp I dp 

This quantity v^^{p) expresses the ratio of the increment of the number of frustrated plaquettes 
to that of antiferromagnetic bonds. 

When p - I, the increment of the number of frustrated plaquettes is larger than that of 
antiferromagnetic bonds. When p =^ 1/2, the converse is realized. There is a turning point 
where the increment of the number of frustrated plaquettes is equal to that of antiferromag- 
netic bonds. It is identified with the point v^'^{p^^) = 1, where 

1 /l\4/3 
/" = 2 + (2) " 0-8969. (6) 

This value is very close to the phase transition point at zero temperature estimated by other 
numerical approaches (Table I). This agreement suggests that the increment of the number 
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Table I. The equations and solutions of v(p) - 1 and the locations pc of the transition points at zero tempera- 
ture in several two-dimensional lattices. The kagome, dual of kagome, extended kagome, and dual of extended 
kagome lattices are depicted in Fig. 3 



Lattice 




Our equation 




Our solution 


Pc 


Square 




2(2/?'*i - 1)3 = 1 




0.8969 


0.8967(1)« 
0.8955(11)" 
0.897(1)1" 
0.8969(1)11 


Triangle 




2(2^*^ - 1)2 = 1 




0.8536 


0.8412(1)^ 
0.833(3)12 


Hexagon 




2(2/5^'='' - 1)5 = 1 




0.9353 


0.9351(2)8 
O.93312 


Kagome 


(2p 


kag _ 1)2 + (2/,k^8 - 1)^ 


= 1 


0.9044 


0.9052(1)8 


Dual of kagome 




2(2;,d-kag _ 1)3 ^ 1 




0.8969 


0.8837(1)8 


Extended kagome 


2(2^ex-kag _ 1)2 + 4(2^ex-kag _ 


1)11=1 


0.9532 


0.9593(2)8 


Dual of extended kagome 




2(2;,d-ex-kag _ 1)2 ^ 1 




0.8536 


0.7948(2)8 



of frustrated plaquettes controls the phases where the system lies. Since temperature has not 
been considered, it is reasonable that the temperature at which our value indicates the transi- 
tion point is zero. Moreover, our argument is consistent with the conjecture on the geometry- 
induced phase transitions at low temperature,^'^ though we have not directly observed the 
entropy of frustration distribution but discussed only the average number of frustrated pla- 
quettes. 

Although our argument yields an extraordinarily accurate correspondence, we cannot re- 
gard it as the derivation of the phase transition point since we have not found any singularity 
in physical quantities. We should examine other lattices for the reasonableness of the agree- 
ment. The cases of other lattices are listed in Table I. The solutions of v{p) = 1 in several 
lattices show fairly good agreement as in the case of the square lattice. These results demon- 
strate that the correspondence between the transition point and the solution of v{p) = 1 is 
not accidental on a particular lattice. In addition, our values in the dual of kagome lattice 
(pd-kag ^ 0.8969) and the dual of extended kagome lattice (^^^-ex-kag ^ 0.8536) correspond 
to the values in the square lattice (p^i = 0.8969) and the triangular lattice (p^ = 0.8536), 
respectively, though the numerical estimates show differences in these lattices. It is noted that 
the dual of kagome lattice and the dual of extended kagome lattice are constructed of squares 
and triangles, respectively, as in Fig. 3. This fact shows that our argument imperfectly reflects 
the global structure of lattices but produces an eff"ective approximation with the local structure 
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Fig. 3. Two-dimensional lattices: (a) The kagome, (b) dual of kagome, (c) extended kagome, and (d) dual of 
extended kagome lattices. 



of lattices. Furthermore, it is remarkable that the result on the dual of extended kagome lattice 
has a larger deviation than the triangular lattice. Our values on the dual of extended kagome 
and triangular lattices are p'^-^^-kag _ ^tn _ 0.8536, whereas the transition points are numeri- 
cally estimated as pf"""''"^ = 0.7948(2)^ and pf = 0.8412(1),^ 0.833(3),'2 respectively. Thus, 
the global structure plays an essential role on the phase transition in the model on the former 
lattice. Indeed, we can find other geometry, for example hexagons, in the dual of extended 
kagome lattice. 

3. Hierarchical lattices 

We next apply the method developed in the last section to the +J Ising model on the 
hierarchical lattices.'^ The lattices are generated by iterating the process that each single 
bond is replaced with a unit of the lattice as depicted in Fig. 4. We can obtain the accurate 
transition point on the hierarchical lattices with renormalization-group calculations.'^''^ The 
lattices are, thus, a good testing ground for the comparison of estimates of the transition 
point. In addition, it is valuable to examine whether or not our argument yields the accurate 
transition points in unusual lattices. We examine the lattices 1, 2, and 3 in Fig. 4. 

Let us calculate the average number of frustrated plaquettes in the hierarchical lattices. 
We first consider the lattice 1. There are plaquettes of various shapes in the lattice, but a 
rule is found. The plaquette of the type generated for the first time by the k-th substitution 
is composed of 2*"^' bonds. For example, a plaquette composed of 8 (= 2^^') bonds appears 
in the second substitution (Fig. 4). The number of plaquettes of 2^"^' bonds in the lattice pro- 
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Fig. 4. The units of hierarchical lattices 1,2, and 3 examined in the present paper and the replacement of each 
single bond with the unit of the lattice 1 . 



duced by iterating the replacement process n times is 4""*", since the number of plaquettes of 
a particular type quadruples in the single replacement. For example, the number of plaquettes 
of 4 (= 2^^') bonds in the lattice after the second substitution is 4 (= 4^'^) (Fig. 4). With 
the periodic boundary condition, we have another plaquette composed of 2""^^ bonds on the 
surface of the lattice. The average number of frustrated plaquettes is, therefore, represented 
as 



^hier. 



\p) = J]4''-W-[l-(2p-ir'}. 



(7) 



/t=i 



The derivative corresponding to Eq. (5) in the infinite- volume limit is 

^^'''-\p) = lim y 2-'^\2p - lf"-\ (8) 

k=\ 

Note that the average number N^^^^'^{p) of antiferromagnetic bonds is 4"(1 - p). Evaluating 
the solution ofv^'^'-^(p^^'-^) = 1 in the case of n = 20, we obtain the value p^'^'-^ = 0.9477. 
Since the same value is obtained even in the case of n = 30, it is reasonable to interpret this 
value as the actual solution. 

In order to compare this value with the location of the transition point at zero temperature, 
we derive the transition point with the renormalization-group method proposed by Nobre. " 
We first prepare a pool of M couplings following the +J distribution, where we set M and 
J equal to 10^ and 1, respectively. A unit of the lattice is built of bonds randomly taken 
from the pool. The unit is renormalized in the ordinary manner of real-space renormalization 
group, which is the reverse operation of the construction of the lattice. Then, a renormalized 
bond is produced and we put it in a new pool. This scheme is executed M times and the 
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Table II. The equations and solutions of v(p) - 1 and the locations pc of the transition points at zero temper- 
ature in the hierarchical lattices in Fig. 4. The indices k of summations in the equations run from 1 to infinity. 

Lattice Our equation Our solution pc 

1 22"*^'(2/''='"'- 1)^'""' = 1 0.9477 0.9215(3) 

2 2 3(1)* (2p- 1)3-2' '-' = 1 0.9089 0.8935(1) 

3 2 3-'^-'{4(2/''='"-3 - l)2-3'-'-i + 8(2/''='"-3 - i)4-3'-'-i) = 1 0.9187 0.8951(3) 



new pool is filled by M renormalized bonds as a result. This is a single renormalization- 
group transformation for the whole system. After the renormalization for the whole system 
30 times, we determine the transition point in terms of the average of resulting couplings. If 
the average is larger (smaller) than 10 (10~^), we regard the system to be in the ferromagnetic 
(paramagnetic or spin-glass) phase. To reduce statistical errors, we run 100 samples. 

The solutions of v{p) = 1 have a bit difference but qualitatively close to the estimates 
of the transition points at zero temperature in the respective lattices listed in Table II. It is 
remarkable that the results correctly express the order of locations of the transition points. In 
particular, a slight difference in the lattice 2 and 3 is distinguished. This result demonstrates 
the effectiveness of our argument to qualitatively predict the transition points also on the 
unusual lattices. 

4. Three-body Ising spin glasses with binary couplings 

In this section we apply our method to a different type of models from previous ones, the 
±7 Ising model with the three-body interactions, 

H = - Yj JijkSiSjSk. (9) 

<i,j,k> 

Couplings Jiji, are governed by the same distribution as the previous models with the two- 
body interactions. One of the reasons to investigate this model is that the order-disorder phase 
boundary in this model is equivalent to the accuracy threshold in a kind of topological quan- 
tum error-correcting codes, namely color codes. ^^' ^^ 

Let us consider the model on the triangular lattice to illustrate the frustration in this model. 
Three spins on the same triangle interact with each other and couplings Jtjk reside on trian- 
gles. A loop in this model is composed of triangles. Each triangle in a loop connects to 
the neighboring triangles on one or two vertices. We call the elementary loop which can be 
frustrated "unit of frustration", which corresponds to a plaquette in the previous cases. The 
unit of frustration in the triangular lattice is a hexagon composed of six triangles depicted as 
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Fig. 5. The unit of frustration in the ±J model with the three-body interactions on the triangular lattice. The 
unit is a hexagon composed of six triangles depicted by thick lines. 



Fig. 5. If the number of negative couplings in a unit of frustration is odd, the unit is frustrated. 
However, this simple rule does not determine whether an arbitrary loop other than a unit of 
frustration is frustrated or not. Even if the number of negative couplings on the loop is odd, 
it is possible that a spin configuration realizes the lower energy state in all couplings. On the 
other hand, whenever there is no frustrated unit in a loop, we can find the spin configuration 
which satisfies all couplings. In other words, loops are frustrated only with frustrated units. 
This fact means that units of frustration are elements to study the nature of frustration in the 
system, and larger loops can miss details of frustration as in the previous models with the 
two-body interactions. 

The average number A^^' '"(p) of frustrated units is represented as 

NyHp) = ^ll-(2p-\)% (10) 

where A'^ denotes the number of spins equal to the number of units of frustration in this lattice. 
The average number N^'^^ of negative triangles is 2N(l - p), and the derivative of N^/^''^{p) 
with respect to A^^' '"(p) is 



V 



v'^''(p) = 3(2p-\f. (11) 

We obtain the value jc^'" = 0.9014, where v^''"(;?^''") = 1. Since the transition point at 
zero temperature in the model, to the best of our knowledge, has not been numerically esti- 
mated in any other studies, we compare this value with the concentration of negative triangles 
at the transition point on the NL instead of the point at zero temperature. The critical con- 
centrations on the NL and at zero temperature are expected to have little difference. Indeed, 
the two values are very close in the case of the Union-Jack lattice, where the transition point 
at zero temperature has been estimated. ^^ The values are compared in Table III. The result 
shows fairly good agreement. Furthermore, it is interesting that our value from v{p) = 1 in the 
triangular lattice is a little smaller than that in the Union-Jack lattice as the critical concentra- 
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Table III. The equations and solutions of v(p) - 1 and the critical concentrations pc and p,„ of negative 
triangles at zero temperature and on the NL, respectively, in the three-body +J model on the triangular and 
Union-Jack lattices. The value pc is expected to be sightly smaller than the value p„ in the triangular lattice as 
in the Union-Jack lattice. 



Lattice 



Our equation 



Our solution 



Pc 



Pm 



Triangle 



3(2/5 



3,tri 



1)^ 



0.9014 



Union- Jack (2p^' ""^ - if + 2(2p^' "-J - 1)7 = i 0.9059 0.8944(1)' 



0.8903(1)20 
0.891(2)1^ 
0.89075(5)20 
0.891(2)21 



tions on the NL in the two lattices which has been derived by the previous study^" accepted 
as a highly accurate prediction. 

Comparing the result for the triangular lattice with the results for the models with the 
two-body interactions on the triangular and hexagonal lattices shown in Table I, we find that 
our argument properly distinguish the number of bodies in the couplings. In particular, the 
difference between the three-body ±7 model on the triangular lattice and the two-body ±7 
model on the hexagonal lattice is worthy of attention, since the units of frustration in these 
two models have the same shape. This coincidence leads to the same term (2p- 1)^ in the both 
functions v(p) accordingly. The coefficients of the terms, 3 in the three-body ±7 model and 2 
in the two-body ±7 model, generate an appropriate difference in the results. The coefficients 
express the increment of frustrated units by an additional negative coupling to the perfect 
pure system. 

5. SK model 

We investigate the standard infinite-range model in spin glasses called the SK model,^ 

H = -Y,JijSiSj, (12) 

where any two spins interact with each other. The couplings 7,^ is governed by the Gaussian 
distribution, 



P(,,)= /i^exp(-^(7,-:^ 



(13) 



where A^ denotes the number of spins, and A^ is introduced in the distribution to make resulting 
extensive quantities be proportional to A'^. 

We consider the number of frustrated plaquettes. A plaquette in this model is a triangle 
composed of three bonds. Since quenched variables 7,y can take arbitrary real values, their 
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product over a plaquette also takes arbitrary values. However, we do not take account of the 
magnitude of the product. Our attention is focused on whether or not a plaquette is frustrated. 
This is one of the ways of generalization of the method developed in the previous cases. The 
average number of frustrated plaquettes is represented as 



Nf(Jo) 



T.\{^-Ul^ 



<ij>ec 



Jo 



(14) 



A^(A^-l)(A^-2) 1^_/^/ ■/( 



6 ^[ \ \-j2N) 

where A^(A'^ - 1)(A'^ - 2)/6 expresses the number of plaquettes equal to the number of combi- 
nations of arbitrary three sites on the lattice, and erf(-) is the error function. 



^'^(ii^r-^-i-^)- 



The average number of antiferromagnetic bonds is 



.f(.) = ^^i(.-e.(4l,. 



In the infinite- volume limit N ^> oo the derivative of Nf^{Jo) with respect to N^^(Jo) is 



V 



/ 

/^^^SK/7„^\-l dl 
,SK 



The solution of the equation v^^(7q^) = 1 is J^^ = ^|n|2. This result does not agree with 
the well known exact transition point at zero temperature Jq = P^ but precisely identical 
with the transition point under the assumption of replica symmetry,^ which is derived from 
the equation of state of the ferromagnetic order parameter m at zero temperature, 

m = erf — — m . (18) 

\V2 / 

This equation has a solution with m i^ if the slope of the function of the right-hand side at 
m = 0, which is ^I/ttJq, is larger than 1 . Therefore, the phase transition point is determined 
by v^^(7o) = 1 that is identical to the condition in our argument. Our result with frustration 
accords with the replica-symmetric solution not only in the solution but also in the equation 
for the determination of the transition point. 

6. Summary and discussion 

We have reported a relationship between frustration and transition point of the Ising spin 
glasses. The relation reveals that the concentration of ferromagnetic bonds in the system is 
rather close to the critical concentration at zero temperature when the derivative of the aver- 
age number of frustrated plaquettes with respect to the average number of antiferromagnetic 
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bonds is equal to unity. This relation is confirmed in various models. In particular, the cases 
in the two-dimensional lattices show good correspondence. In the SK model, the solution 
derived from our argument exactly agrees with the transition point under the assumption of 
replica symmetry, and the slope of the average number of frustrated plaquettes corresponds 
to the slope of the function in the equation of state of the ferromagnetic order parameter at 
zero temperature. 

Although we have detected no conventional sign of phase transitions, the agreements 
between our values and the transition points are naturally regarded not to be accidents. One 
of the reasons is that our argument attains pretty accurate predictions in most of the models 
which we have applied the method to. The variety of models which the agreements are found 
is especially remarkable. Our argument is, thus, expected to yield approximate locations of 
the transition points. 

Our approach to spin glasses by means of frustration is novel and different from con- 
ventional ones. The representative conventional study is an attempt to characterize phases 
in terms of the change in the free energy by increase of frustration.^^ Both our study and the 
conventional one focus on the change of frustration to discuss phases. Hence, they might have 
a deep connection, but we have not found it yet. It is notable that our approach enables us 
to extract the approximate locations of the transition points specifically, which is difficult for 
conventional ones. 

Furthermore, it is remarkable that the result for the SK model, which corresponds not to 
the exact solution but to the replica symmetric one, shows deep significance of the deviation 
from the exact solution. This result suggests that we should not dispose the gap between 
our solution of v = 1 and the exact one as a simple defect of our method. Moreover, the 
correspondence of v^^(7o) and the slope of the equation of state provides an insight into our 
heuristic condition v = 1. Also in other cases, v is expected to be concerned with the slope of 
the equation of state under some assumption like that of replica symmetry. 

Simplicity of our calculation in contrast with the accuracy of the results has practical util- 
ity in the context of topological quantum error-correcting codes, where the accuracy threshold 
corresponds to the phase boundary of spin glasses. Our method is useful for the first approxi- 
mation of the accuracy threshold. The simplicity, moreover, illuminates the unknown elegant 
physics in spin glasses, which is regarded as a complex subject usually. 

Since we have considered only frustration to discuss phase transitions, our result seems to 
support the conjecture on the geometry-induced phase transition.^' ^ Nevertheless, our values 
are different from expected locations from the conjecture. ^'^'^^"^^'^"^"^^ The reason is that we 
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have observed the number of frustrated plaquettes, whereas the entropy of frustration dis- 
tribution has been considered in the conjecture.^' ^ Our result suggests the existence of the 
geometry-induced phase transition, but it does not support the proposition that the root of the 
transition at zero temperature is a singularity in the entropy of frustration distribution. 

Since our results have the suggestion about the relationship among frustration, phases, 
and the replica symmetry, the origin of the agreements might lead to the fruitful structure 
behind spin glasses in general. The origin is, however, not obvious yet. To clarify it, further 
investigation is necessary. In particular, the lack of the effect from the global structure of 
lattices and from the nature of spins should be resolved. It is especially important to clarify 
the reason why our results are consistent with the transition points of models with the Ising 
spins. 
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